THE NUMBER OF RHOMBUS TILINGS OF A SYMMETRIC 
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SYMMETRY AXIS, II 
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Abstract. We compute the number of rhombus tihngs of a hexagon with side lengths 
N, M, N, N, M, N, with N and M having the same parity, which contain a particular 
rhombus next to the center of the hexagon. The special case iV = M of one of our 
results solves a problem posed by Propp. In the proofs, Hankel determinants featuring 
Bernoulli numbers play an important role. 



1. Introduction 



Let a, b and c be positive integers, and consider a hexagon with side lengths a, b, c, a, b, c 
whose angles are 120° (see Figure |l].a). The subject of enumerating rhombus tilings of this 
hexagon (cf. Figure |I[b; here, and in the sequel, by a rhombus we always mean a rhombus 
with side lengths 1 and angles of 60° and 120°) gained a lot of interest recently. This 
interest comes from two facts. First, it is a rich source of non-trivial enumeration problems 
which have (or appear to have) beautiful solutions (see e.g. |^, ||, ^ 0, |2^, 0). 

Second, these problems are very often related to the theory of symmetric functions and/or 
the representation theory of classical and quantum Lie algebras, and to statistical physics 
(sometimes in disguise; see e.g. H, 0, |2|, |2|, ||, ||, |9Hg]). 

As is well-known, the total number of rhombus tilings of a hexagon with side lengths 
a, b, c, a, b, c equals 



i=i j=i k=i •' 



:i.ii 



(This follows from MacMahon's enumeration Sec. 429, g ^ 1; proof in Sec. 494] of 
all plane partitions contained in an a x 6 x c box, as these are in bijection with rhombus 
tilings of a hexagon with side lengths a, 6, c, a, 6, c, as explained e.g. in [p]] .) 

A natural question to be asked is what the distribution of the rhombi in a random 
tiling is. On an asymptotic level, this question was answered by Cohn, Larsen and Propp 
TOl) . On the exact (enumerative) level, the significant contributions are 0, |1^, |16|. The 



most general result was obtained in by the authors, where the number of all rhombus 
tilings of a hexagon with side lengths N, M, N, N, M, N was computed, which contain an 
arbitrary fixed rhombus on the symmetry axis which cuts through the sides of length M. 
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The purpose of this paper is to add other results in this direction. We compute the 
number of all rhombus tilings of a hexagon which contain a particular rhombus on the 
"other" symmetry axis, i.e., the symmetry axis which runs in parallel to the sides of length 
M. In difference to we are not able to solve this problem for an arbitrary rhombus 
on this symmetry axis, but only for rhombi which are close to the center. In fact, the case 



of the central rhombus is already covered by the papers m O, 16]. We provide results for 



the "next" three cases, i.e., for a rhombus which is by one, two, or three "half units" off 
the center, see Theorems below. In fact, as our proofs show, the computations become 
increasingly harder, and more elaborate, as we move the rhombus gradually farther away 
from the center, so that it seems highly unlikely that a uniform formula, similar to the 
one in for the other symmetry axis, can be found for an arbitrary rhombus. (For 



further comments on this issue see Section ||, (2).) 

Here are our results. (Each of the following formulas has to be interpreted as the 
appropriate limit if singularities are encountered. For example, if we directly set m = 1 
in ( |1.6| ) then the term (m — 1) in the second line becomes 0, and, on the other hand, we 
have a singularity in the sum in the third line caused by the term (3 — n — m)/j in the 
summand for h = n — 1. The correct way to interpret the expression is as the limit as m 
goes to 1.) 



Theorem 1. Let n and m be positive integers. The number of rhombus tilings of a 
hexagon with side lengths 2?7,, 2m, 2n, 2n, 2m, 2n, which contain the rhombus above and 
next to the center of the hexagon [see Figure^, a; the rhombus which is contained in every 
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a. A hexagon with sides N, M, N, 
N, M, N, where N = 4, M = 2. 



b. A hexagon with sides A^, M, A^, 
A^, M, N, where A^ = 3, M = 3. 



Figure 2. 



tiling is shaded), equals 

\ n / \ n / \ m / 



(4n+2m-l\ 
\ 2n+m ) 



+ 



An + 2 



{n + m)^ {n + l)(2r?, — l)(n + m — l)(n + m + 1) 



X 



>^ {2)h (1 - n)h (I + n)/i (1 - n - m);, (1 + n + m)h 
^ (1);, (2 + n)h (f - n)42 + n + m)/, (2 - n - m)h 



2n 2m 2n 



i+j+k-l 



(1.2) 

where the shifted factorial (a)^ is defined by {a)k := a(a+l) ■ ■ ■ (a+A; — 1), k > 1, (a)o := 1. 



xnnn,+,+,_2- 

1=1 j=i fc=i 



Theorem 2. Let n be a nonnegative integer and m be a positive integer. The number of 
rhombus tilings of a hexagon with side lengths 2n + l, 2m — 1, 2n + l,2n + l, 2m — 1, 2n + l, 
which contain the rhombus above and next to the center of the hexagon {see Figure ||. 6; 
the rhombus which is contained in every tiling is shaded), equals 



Mn+2m\ 
V 2n+m / 



+ 



An 



{n + m)^ {n + l){2n — l){n + m — l){n + m + 1) 



n-l 



X 



>p (2)/, (1 - n)h (f + n)h{l - n - m)h {I + n + m)h 
^ {l)h (2 + n)h (I -n)h{2 + n + m)h {2 - n - m)h 



h=0 



2n+l 2m- 1 2n+l 

n n n 

i=i j=i k=i 



i+j+k-l 
i+j+k-2 



■ (1-3) 
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Theorem 3. Let n and m be positive integers. The number of rhombus tilings of a 
hexagon with side lengths 2n, 2m — 1, 2n, 2n, 2m — 1, 2n, which contain the rhombus above 
and next to the central rhombus (see Figure^. a; the rhombus which is contained in every 
tiling is shaded), equals 

(2--i)(r_t)(")e::i) 



(n + m - 1) (n + m) 

n {n + 1) {2n — 3) {2n — 1) (m^ — m + n + 2m n + + 1) 
{n — 1) {n + m — 1) (ra + m) {2n + 1) 

6 (3)/. {\)h (2 - n)n (I + n\{2-n- m)^ (1 + n + m\ 



(n + m - 2)(n + m + 1) ^ (1)^ (2 + n)h (| - n)/, (2 + n + m)h (3 - n - m\, 

2n 2m~l 2n 



i=i j=i k=i 

Theorem 4. Let n and m be positive integers. The number of rhombus tilings of a 
hexagon with side lengths 2?t, — 1 , 2m, 2n — 1 , 2n — 1 , 2m, 2n — l, which contain the rhombus 
above and next to the central rhombus (see Figure^.b; the rhombus which is contained in 
every tiling is shaded), equals 

(n + m - 1) (n + m) 

n {n + l) {2n - 3) (2n - 1) (m^ - m - 3n + 2mn + + 2) 
(n — 1) (n + m — 1) (n + m) (2n + 1) 



+ 



(3)^ (I);, (2 - n\ (I + n)h{2-n- m)^ (1 + n + m)^ 



(n + m - 2)(n + m + 1) ^ (1)^ (2 + n)^ (| - n)42 + n + m);, (3 - n - m)^, 

2n-l 2m 2n-l 



^ n n n <^-^' 

i=i j=i k=i 

Theorem 5. Let n and m be positive integers. The number of rhombus tilings of a 
hexagon with side lengths 2n, 2m, 2n, 2n, 2m, 2n, which contain the rhombus above and 
next to the rhombus which is adjacent to the center of the hexagon (see Figure 0.a; the 
rhombus which is contained in every tiling is shaded), equals 

/n+m-l\ f2n+2\ /2n+m-l\ /2n+ni\ 
V m Ara-lA n / V 2n+l / 

2{2n - 3){2n - l){2n + 2){n + m - l){n + m + ('"+'r"') 
2 (n + 2) (n + 3) {2n - 1) {2n - 3) X(m, n) 24(m - l)(2n + m + l)(2n + l)(2n + 3) 
(ra + m - 1) (n + mf (n + m + l) {n + m - 2){n + m + 2) 

^ (4),(l-n),(| + n),(2-n-m),(2 + n + m), \ ^^^ ^+J + k-l 
^(l),(4 + n),(|-n)43 + n + m)43-n-m),;l|ljllz+j + A;-2' 



X 
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a. A hexagon with sides A^, M, N, b. A hexagon with sides A^, M, N, 

N, M, N, where A^ = 2, M = 3. N,M, N, where A^ = 3, M = 2. 

Figure 3. 

where 

X{m,n) = — l+4m —3m —n+8mn+llm n — 12171 n — 6m n+8n +22mn —4m n — 24m n 
- 2m^n^ + 15n^ + 16mn^ - 28m^n^ - 8m^n^ - 5n^ - 8mn^ - 12m^n'^ - Un^ - 8mn^ - 2n^. 



Theorem 6. Let n and m he positive integers. The number of rhombus tilings of a 
hexagon with side lengths 2n + 1, 2m — 1, 2n + 1, 2r;, + 1, 2m — 1, 2?7, + 1, which contain the 
rhombus above and next to the rhombus which is adjacent to the center of the hexagon 
(see Figure ^.b; the rhombus which is contained in every tiling is shaded), equals 

\ m }\ m-1 ) V n-1 / V n+1 / 

(m + l)(m + 2)(2n - 3)(n + m - + m + 1) (^t^) (f+D 

/ (n + 2) (n + 3) (2n - 1) (2n - 3) Y{m, n) 24(m - l)n{2n + m + l)(2n + 3) 
y {n + m-l){n + mf (n + m + l) {n + m - 2){n + m + 2) 

(4), (1 - n), (I + n), (2 - n - m), (2 + n + m), \ t + j + k-1 

■^(l),(4 + n),(|-n),(3 + n + m),(3-n-m),; A| ii^ + j + /c-2' 

(1.7) 

where 

y (m, n) = 1 — 2m + m + 5n — 4mn — 2m n + 4m n — 3m n + 22n — 4mn — 4m n 
-l2m^n^-2m'^n^+bQn^-lQmn^-2Qm^n^-8m^n^+2,^n'^-28mn^-l2m'^n^+bn^-8mn^-2n^. 

In general, the sum in (|1.2|) and (|1.3|) (note that it is indeed exactly the same sum), 
the sum in ( |1.4| ) and ( |1.5| ) (it is indeed exactly the same sum), and the sum in ( |1.6| ) and 
( [1.7| ) (again, it is indeed exactly the same sum), does not simplify. However, in the case 
that n and m are roughly of the same size, the sum does simplify. For the sake of brevity, 
we give here just a sample of corollaries to Theorems |1| and 0, the first two statements of 
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Figure 4. 




a. A hexagon with sides A^, M, A'", 
A^, M, N, where N = 4, M = 2. 



M M 




b. A hexagon with sides N, M, N, 
N, M, N, where A^ = 3, M = 3. 



which solve the second part of Problem 4 in Propp's list [0. We wish to emphasize that 
there are similar corollaries to Theorems 



Corollary 7. Let n be a nonnegative integer. The number of rhombus tilings of a hexagon 
with all sides of length 2n, which contain the rhombus above and next to the center of the 
hexagon, equals 

f2n\3\ 2n 2n 2n . . J ^ 



1 (17 



12 C) 



i+j+k-2 



i=i j=i k=i 

The number of rhombus tilings of a hexagon with all sides of length 2n+ 1, which contain 
the rhombus above and next to the center of the hexagon, equals 

('2r^^3 \ 2n+l 2n+l 2n+l ■ , ■ , ; 

n n n Sj^l^- 

i=i j=i k=i •' 

The number of rhombus tilings of a hexagon with side lengths 2n, 2n + 2, 2?2, 2n, 2?t, + 2, 2n, 
which contain the rhombus above and next to the center of the hexagon, equals 

(2n\^ \ 2n 2n+2 2n 




l_(10n + 2)_gT 

3 (6- +3) e:D, 



The number of rhombus tilings of a hexagon with side lengths 2n + l,2n — l,2n + l,2n + 
1, 2n — l,2n + 1, which contain the rhombus above and next to the center of the hexagon, 
equals 

1 I Q\ ('2n\3\ 2n+l 2n-l 2n+l • , ■ , , 

Mi5!i±^4i^ n n n ^^j^l^- a-") 



24n 



\3n) 



i=l j=l k 

The number of rhombus tilings of a hexagon with side lengths 2n + 2, 2n, 2n + 2, 2n + 
2, 2n, 2n + 2, which contain the rhombus above and next to the center of the hexagon. 
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equals 

2n+2 2n 2n+2 ^ 

n n n Si^^- f^-"' 

T/ie number of rhombus tilings of a hexagon with side lengths 2n + 3,2n — 1, 2n + 3, 2ri + 
3, 2n — l,2n + 3, which contain the rhombus above and next to the center of the hexagon, 
equals 

<^(r 2 , n , o\ /an^S \ 2n+3 271-1 2n+3 • , ■ , , 

2(6n^ + 9n + 2) (J \ + 

(^ + 1)2 ^6n+4N I 11 11 11 , + ■ + fc_2- ^'•'''^ 

^ ' \Zn^2l / i=\ j=\ k=l 

If m is not ra + l, n or n — 1, then one may still try to obtain at least estimates 
for the number of rhombus tilings that contain this particular rhombus. Indeed, from 
Theorems we are able to derive an "arcsine law" for this kind of enumeration, which 
is analogous to the ones in 0] and |]T3[ . 




Corollary 8. Let a be any nonnegative real number, and consider a hexagon with side 
lengths 2n, 2m, 2n, 2n, 2m, 2n. For m ~ an, the proportion of the rhombus tilings that 
contain the rhombus above and next to the center of the hexagon in the total number of 
rhombus tilings is ~ arcsin(l/(a + 1)) as n tends to infinity. The same is true for 
a hexagon with side lengths 2n + 1, 2m — l,2n + l,2n + 1, 2m — 1, 2n + 1. Moreover, 
for hexagons with side lengths 2n, 2m — 1, 2n, 2n, 2m — l,2n or hexagons with side lengths 
2n—l, 2m, 2n — l,2n—l, 2m, 2n—l, the same result holds for the proportion of the rhombus 
tilings that contain the rhombus above and next to the central rhombus in the total number 
of rhombus tilings, as well as for the proportion of the rhombus tilings that contain the 
rhombus above and next to the rhombus which is adjacent to the center of the hexagon in 
the total number of rhombus tilings of hexagons with side lengths 2n, 2m, 2n, 2n, 2m, 2n 
or hexagons with side lengths 2n + 1, 2m — l,2n + 1, 2n + 1, 2m — 1, 2n + 1. 

Also this result is (as well as Corollary 4 in and Theorem 1.3 in [|13|; see the respective 



comments in [Q) in accordance with Conjecture 1 in |T0[, to which it adds evidence in 
further special instances. 

In the next sections we describe proofs of Theorems and of Corollaries |^ and 
^. In Section 2 we provide proofs of Corollaries |^ and ^ and we outline the proofs of 
Theorems ^|-|^, the latter consisting of two basic steps. In the first step we build on 
the approach of Helfgott and Gessel in [|l^ , a short summary of which is the contents of 
Section 3. It allows to write the number that we are interested in in form of a determinant. 
The evaluation of this determinant is not easy and is carried out in detail in Section 4. 
For the evaluation we follow a "method" that was first introduced in p2| (see the tutorial 



description in p3|. Sec. 2.4] or [0, Sec. 2]). For accomplishing the required computations, 
we need to evaluate certain Hankel determinants featuring Bernoulli numbers, which are, 
in fact, of independent interest. As it turns out, some of the evaluations of these Hankel 
determinants are already known, provided certain results about orthogonal polynomials, 
in particular, about continuous Hahn polynomials, and continued fractions are properly 
combined. For the convenience of the reader, we collect these facts, and their implications, 
in Section 5. In particular, the evaluation of the relevant Hankel determinants featuring 
Bernoulli numbers is given in Theorem ^ However, in the proofs of Theorems § and ^ 



(more precisely, in the proof of the subordinate Lemma O) we encounter a certain Hankel 
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determinant of Bernoulli numbers (see ( |6.1| )), the evaluation of which requires considerable 
effort. (This is one of the added difficulties mentioned earlier in comparison to the proofs 
of Theorems |I| and |^.) We evaluate this Hankel determinant by combining the knowledge 
about continuous Hahn polynomials with a recent theorem on orthogonal polynomials due 
to Leclerc (restated here as Theorem and applying some integral calculus (see the 
proof of Lemma ^Bp. Section 6 is devoted to provide the details of these calculations. In 
Section ^ we make explicit a few unusual evaluations of Hankel determinants of Bernoulli 
polynomials, which are implicit in the proofs of our enumeration results. Finally, in 
Section ||, we point to further directions in this research, and propose a few open problems. 

2. Outline of proofs 

Here we outline the proofs of Theorems and we deduce Corollaries ^ and ||. We 
fill in the details in the subsequent sections. 

Proof of Theorems Following the approach of Helfgott and Gessel Q (see 
Section 3), we may write the number of rhombus tilings of a hexagon with side lengths 
N, M, N, N, M, N, which contain an arbitrary rhombus on the [N + M)-long vertical 
symmetry axis (see Figure ^), in form of a determinant. This determinant is given by 
Proposition |TT]. That is to say, in order to prove Theorems ^ and |^, we need to evaluate 
the determinant in (|3l^ ) with / = , and in order to prove Theorems || and |, we need 
to evaluate the same determinant with / = i^±M±l^ and in order to prove Theorems ^ 
and ^ we need to evaluate the same determinant with / = K±M+^_ Modulo replacement 
of parameters, we may thus concentrate on the determinants D{n,n, N), D{n,n — 1,N) 
and D{n,n — 2, N), where 

D{a,b,N):=Jet^('^M. (2.1) 

^ s=—a ' 

The evaluation of determinants -D(n, n, A^), D{n,n — 1, N) and D{n,n — 2, N) is then 
carried out in Lemmas [12], and |T^, respectively. 

Theorem Q follows upon combining Proposition 11 with N = 2n, M = 2m, I = n + m 
and Lemma O with n replaced by n + m and N = 2n — 1, and some rearrangement of 



terms. Similarly, Theorem y follows upon combining Proposition [Tl| with = 2r;, + 1, 
M = 2m — 1, I = n + m and Lemma with n replaced by n + m and = 2n. 

Likewise, Theorem follows upon combining Proposition O with N = 2n, M = 2m — 1, 



I = n + m and ([4.18|) with n replaced by n + m and m replaced by n, and Theorem ^ 



follows upon combining Proposition 11 with N = 2n — 1, M = 2m, I = n + m and (^4.19|) 
with n replaced by n + m and m replaced by n — 1. 

Likewise, Theorem ^ follows upon combining Proposition |Tl] with = 2n, M = 2m, 
I = n + m + 1 and ( ^4.23|) with n replaced by n + m + 1 and m replaced by n — 1, and 
Theorem ^ follows upon combining Proposition |ll] with A^ = 2n + 1, M = 2m — 1, 



I = n + m + 1 and ( [4.24| ) with n replaced by n + m + 1 and m replaced by n. □ 



Proof of Corollary 0. We have to compute the value of the expressions ( pr2D and 
]3]) for m = n + 1 (in order to establish ( |1.10| ) and ( |1.9| )), m = n (in order to establish 
^ and ( p.. lip ), and m = n — 1 (in order to establish ( |1.12| ) and ( |1.13|) ). Clearly, except 



for trivial manipulations, we will be done once we are able to evaluate the sum in (|T 
and (|1.3|) (it is indeed the same sum!) for m = n + 1, m = n, respectively m = n — 1. 
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We treat the case m = n first. We claim that 

n— 1 



. _ fi 



(2);, (3/2 + n)h (1 - n)h (1 + 2n)h (1 - 2n),, _ (n + l)(2n - if ( 1 1 © 



+ 



'(2.2) 

Let us denote the sum by S{n) and its summand by -F(n, i). We use the Gosper-Zeilberger 



h=0 



algorithm [5^, |51| to obtain the relation 

Qn^{n + 2){<6n + l)(6ra + 5) F{n, i) - 6(n + l){2n - lf{?,n + l)(3r2 + 2) F{n + 1, i) 

= G{n,i + l)-G{n,i), (2.3) 

with 

nin + 2)i2n - 2h - l){2n - h ~ I) 
^ ' ^ (/i + l)(/i-n)(2n-/i)(2n-/i + l)(2n + /i + 2) 
X (144n^ - 432/i^n^ - 432/in^ + 312n^ - 936/i^n^ - 936/in^ + 236n^ + 

+ 216/1^^2 - 588/1^^2 - mQhv? + TOn^ + mh^n + 234/i^n - 83/i^n - 200/in 

+ 6n + 24/i^ + 48/i^ + G/i^ - 18/i)F(n, i). 

Summation of the relation (|2.3|) from z = Otoi = ra — 1, and little rearrangement, leads 
to the recurrence 

Qn^{n + 2){<on + l)(6n + 5) S{n) - 6(n + l)(2n - l)2(3n + l)(3n + 2) S(n + 1) 

_{n + 2){2n - 1)2(36^3 + eOn^ + 29n + 3) 

~ 2{n + l) ■ 

for the sum in (|2.2| ). (Paule and Schorn's Mathematica implementation of the Gosper- 
Zeilberger algorithm, which is the one we used, gives this recurrence directly.) Since 
S{1) = 1, and since the right-hand side of (|2.2|) satisfies the same recurrence, equation 
( p.2|) is proved. 

The procedure in the other two cases is analogous. The respective evaluations that we 
need to prove are 

(2)fe (3/2 + n)H (1 - n)h (2 + 2n)fe {-2n)H ^ {n + l)^{2n-l) ( _2 I QZ'^ 

h=0 



il)U3/2-n)f,{2 + n)Ul-2n)f,{3 + 2n)h (2n + 1)^ ^3 3 (2")^ 



:2.4) 



in the case that m = n + 1, and 



n-l 



(2);, (3/2 + n)^ (1 - 72)42^)/, (2 - 2n)„ 2n(n+l)/ , n £J) 



1 



^^^il)hi3/2-n)h{2 + n)h{3-2n)hil + 2n)H 2n + 1 ^ 12(2r2 - 1) (2n-2^j3 

^ (2.5) 

in the case that m = n — 1. We leave it to the reader to fill in the details. □ 

Proof of Corollary |8|. We concentrate first on the case of a hexagon with side 
lengths 2n, 2m, 2n, 2n, 2m, 2n and the rhombus above and next to the center of the 
hexagon. 

From MacMahon's formula (|1 . 1| ) for the total number of rhombus tilings together with 
Theorem |1| we infer that the proportion of the rhombus tilings that contain the rhombus 
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above and next to the center of the hexagon with sides lengths 2n, 2m, 2n, 2n, 2m, 2n in 
the total number of rhombus tilings is given by 



nm 



f2n\ f2n-l\ f2m\ 
\ n / \ n / \ m / 



1 



f4n+2m-l\ 
V 2n+m ) 



+ 



An + 2 



{n + m)^ {n + \){2n — l)(n + m — l)(n + m + 1) 



X 



n-l 



{2)h (3/2 + n);, (1 -n)h{l + n + m)h {I - n - m)h 



h=0 



{l)h (3/2 - n)h (2 + n)h{2-n- m)h i2 + n + m)h ' 



Using the standard hypergeometric notation 



s 



ai, 
bi, 



; z 



A:=0 



k\ {h 



0'r)k k 



Ik 



iP. 



's k 



(2.6) 



we may write the above expression as 



nm 



1 



Mn+2m-l\ 
V 2n+m I 



+ 



An + 2 



{n + m)^ {n + \)(2n — l)(ra + m — l)(n + m + 1) 
2, 2, 1 + n + m, 1 — — m, 1, 1 + n, 1 — n 



X 



3- ;l 



1,2 — n — m, 2 + r;, + m, 2,| — n,2 + r;, 



[2.1] 



Now we observe that the sum in this expression is in fact a very-well-poised hypergeo- 
metric series, and therefore we can transform it into a more convenient form by applying 
Whipple's transformation (see (2.4.1.1)]) 



7^ 6 



a, 1 + |, 6, c, (i, e, — 
,|,l + a — 6, 1 + a — c, 1 + a — (i, 1 + a — e, 1 + a + A^' 



(a + l)Ar (a — d — e + 1^ 



TV 



a — d -|- l)Ar (a — e + 1' 



4-r 3 



a — 6 — c + 1, c?, e, —N 
a — 6 + l,a — c+1, —a + + e — 



;1 



:2.8) 



(where is a nonnegative integer) with a = 2, 6=l + n + m, c=l — n — m, (i = l, 
e = 3/2 + n, and = n — 1 to it. Thus we obtain the expression 



/4r),+2m-l\ 
V 2n+m ) 



1 



+ 



2n + l 



{n + m)^ {n + m — l){n + m + \) 



X 4^3 



1, 1, 1 + ra, 1 - n 
2 — 72 — m, 2 + n + m, I' 



:2.9) 



Now we substitute m an and perform the limit n ^ oo. The asymptotics of the 
binomials appearing in front of the expression (|2.9| ) is easily determined by means of 
Stirling's formula. For the 4F3-series itself, we may exchange limit and summation by 
uniform convergence. 



lim 4F3 



l,l,|+n,l 



n 



2 — n — an, 2 + n + an, | 



3,1 



1 



1,1 

. 1 '{a + iy 



(2.10) 
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Combining all this, and making use of the identity (see [|^, p. 463, (133)]) 

arcsin \fz 



3 1 Z 
. 2 



(2.11) 



^z{l - z) 

we obtain exactly ^ arcsin(l/(a + 1)) as the asymptotic form of ( [^.91 ), and, hence, as the 
asymptotic form of the proportion of rhombus tilings in the statement of the corollary, as 
desired. 

The case of a hexagon with side lengths 2n + \, 2m — l,2n + l,2n + \, 2m — 1, 2n + 1 
can be handled in (almost) the same way because the sums in (|1.2| ) and ( p..3|) are exactly 
the same. 

For the next two cases, i.e., in order to estimate (|1.4|) and (|1.5|) , we proceed in a similar 
way. Again we write the sum which appears in ( |1.4|) and (|1.5| ) as a yFg-series, 

'3, 1, 1 + n + m, 2 - n - m, 2, 1 + n, 2 

m,2,| 



7-^6 



5 

^3 



n 



m,2 + n 



n 



;1 



2," I 'f I "'';^)2 n^2 -\- n 

apply Whipple's transformation ( |2.8| ), and then let n tend to infinity. Here, the 2-^1-series 
which is obtained is a slightly different one as before (compare ( |2.10| )), 

1 



1,2 

5 , 

L 2 



(a + 1)2 

In order to be able to use (|2.11|) , we use the relation 



1,2 1 



3(a + 1)^ 



2F1 



1,1 
. 2 



(a + i; 



The computation is then completed by straightforward use of Stirling's formula, and 
subsequent simplification. 

The remaining two cases, i.e., the estimations of ( |1.6| ) and ( |1.7| ), can be dealt with in 
just the same manner. The 2-Fi-series which is obtained here is again slightly different. It 
reads 



2F1 



1,2. 

3 , 
. 2 



a + 1 



r 



iF^ 



1,1. 

1 , 

L 2 



2 2 

To the 2-^1-series on the right-hand side we apply the formula (see 

1 



I p. 464, (138)]) 



2Fl 



1 , z 

L 2 



+ 



z arcsm ^z 



Again, the computation is then completed by straightforward use of Stirling's formula, 
and subsequent simplification. □ 

3. From rhombus tilings to determinants 

This section is entirely based on ideas by Helfgott and Gessel |jl6|. These allow us to 
find a determinantal expression for the number of rhombus tilings of a hexagon with side 
lengths N, M, N, N, M, N that contain an arbitrary fixed rhombus on the (A^ + M)-long 
symmetry axis (the vertical symmetry axis in Figure |^). We shall state two auxiliary 
results (Propositions |^ and |10]) without proof (the reader can find the details in |1^), and 
then derive Helfgott and Gessel's determinant (Proposition [T^). It is the specialization 
/ = (A^ + M)/2 of Proposition |ll] (compare the paragraph above (|2.1| )) which in the long 
run leads to a proof of our Theorems ^ and |^, it is the specialization I = [N + M + l)/2 
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Figure 5. 




Figure 6. 



of Proposition 11 which in the long run leads to a proof of our Theorems ^ and ^, and it 
is the specialization / = {N + M + 2)/2 of Proposition |11| which in the long run leads to 
a proof of our Theorems ^ and ^. We do want to alert the reader that we use a different 
convention in our figures of how to draw the hexagons than Helfgott and Gessel. To be 
precise, our figures turn into those in ||T6[ by a rotation by 90°. 

The first observation is that for any rhombus tiling of the hexagon with side lengths 
N, M, N, N, M, N, there are exactly rhombi of the tiling that are cut in two by the (A^+ 
M)-long symmetry axis. Removing these rhombi, and cutting the hexagon in two along 
the symmetry axis, leaves two symmetric halves of trapezoidal shape with "dents". 
The following statement counts rhombus tilings of such "dented trapezoids". (This is 
Lemma 2 in 111 
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Figure 7. 



Proposition 9. The number of rhombus tilings of a semi-hexagon with side lengths N, 
M,N (i.e., the "half" of a hexagon with side lengths N,M,N,N,M,N) and N "dents" 
at positions < Tq < ■ ■ ■ < Ttv-i < N + M (see Figure where N = 4, M = 3, and the 
"dents" are at positions 0, 1, 4 and 5) is 

\i=l / 0<j<j<Af-l \i=l / ~ 

where pi{x) is an arbitrary monic polynomial of degree i in x. (The "standard" case would 
be pi{x) = x*, which corresponds to the Vandermonde determinant.) 

From this proposition, Helfgott and Gessel deduce another enumeration result. (This 
is Proposition 4 in ||16|| .) 

Proposition 10. Let L be a subset of 0,1, N + M — 1 of cardinality at least N . Then 
the number of rhombus tilings of a hexagon with side lengths N, M, N, N, M, N , in which 
the set of rhombi on the {N + M)-long symmetry axis is a subset of L (given in terms 
of the numbers of the rhombi, where the rhombi on the symmetry axis are numbered from 
bottom to top asO, 1, . . .,N+M-1; see Figure^, where N = 4, M = 3, L = {0, 1,3,4,5}; 
the set L consists of the shaded rhombi), is 

(n-ij)„,d|Vi(g''-wa(^')' 

where, again, pi{x) is an arbitrary monic polynomial of degree i in x. 

From this proposition we can derive the following determinantal expression for the 
number of rhombus tilings of a hexagon with side lengths N, M, N, N, M, N which contain 
an arbitrary fixed rhombus on the (A^ + M)-long symmetry axis (see [|T^). 

Proposition 11. The number of rhombus tilings of a hexagon with side lengths N, M, N, 
N,M,N which contain the l-th rhombus, 0<1<N + M-1, on the {N + M)-long 
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symmetry axis (see Figure ^ the rhombus which is contained in every tiling is shaded, 
i.e., 1 = 2), is 

Proof. Let us first count the complementary set, i.e., the rhombus tilings which do 
not contain rhombus /. Obviously, we get this number from Proposition |10| with L = 
{0, 1, . . . , + M — 1} \ {/}. Then, if in addition we choose Pi{x) = {x — ly, by formula 
this number is 

We have to subtract this number from the total number of possible rhombus tilings. This 
can be again expressed by making use of Proposition |10|, this time with L = {0, 1, . . . , + 
M — 1}. So, choosing Pi{x) = (x — ly in (p.l|) again, we obtain 

n^)„,-.(E(^-o-)^(n^).,-^( £ (3.4) 

for the total number of rhombus tilings. 



It should be observed that the determinants in ( |3.3D and ( |3.4| ) are almost the same. 
The only difference is the (0, 0)-entry, which is by 1 less in ( |3.3|) than in ( |3.4| ). Therefore, 
if we expand the determinant in (p.3|) with respect to the first row, we may rewrite (|3.3|) 

as 

The number of rhombus tilings that we are interested in is the difference of (|3.4]) and 
(PD, which is exactly (pj). □ 

4. Determinant evaluations 

Lemma 12. Let n and N be positive integers. Then the determinant D{n,n, N) , as 
defined in ( |2.1| ), is equal to 

[N/2i 



Y[ {{n' - tY~''^\n' - (z - 1/2)2)^-2.+!^ 

2^^iV! {N + iy{n- \N/2] )21n/2]+i ITU 



n 

i=l 

X 



n[N/2\\^ \N/2V.'UZ''' 



, 4(iV + 2)n^ [N/2\ \N/2-\ 

' ' ^ {N -l)N{n^-l){[N/2\+l){\N/2-\+l) 

^""^^ (2), (1 - \N|2^), (3/2 + riV/21 ), (1 - n), (1 + n), . 
^ (1), (2 + [iV/21 ), (3/2 - riV/21 ), (2 + n\{2-n)J- ^ ' > 
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Proof. We proceed in several steps. An outline is as follows. In the first step we make 
the obvious observation that D{n,n,N) is actually a polynomial in n, of degree at most 
A^(A^ + 2). Next, we show that D{n, n, N), as polynomial in n, has a lot of linear factors. 
More precisely, in the second step, we show that is a factor of D{n,n, N). Then, in 
the third step, we show that n«-^i^^(^^ ~ ^2^Af-2i+i jg factor of D(n,n, N). Moreover, 
in the fourth step, we show that nl=i^^('^^ ~ ~ 1/2)^)^"^*"'"^ is a factor of D{n,n, N). 
From a combination of these four steps we are forced to conclude that 

D{n, n, iV) = ( J] {n' - zY-''+\n^ - (z - l/2)Y-^'+' J ■ P{n, N), (4.2) 
^ i=i ^ 

where P{n,N) is a polynomial in n of degree at most 2 \N/2]. Finally, in the fifth step, 
we evaluate P(n, A^) at n = - [A^/2] , - [A^/2] + 1, . . . , [A^/2] - 1, [A^/2] . Namely, we 
show that 



^ ^7 + 1)^ i^ + 2) 



N-i 



]-[W2j ^2{N-2',+l) _ 1/2)2(^-2^+1) ii l(2z + 1) (2z + 2)^ {2l + 3) 



(4.3) 



and, for 1 < e < [A^/2] , that 



P{±e,N) 



^_^YN-2e+l){N~2e)/2 2^^-26+1 Y[ 



2e-l^,2 



L7V/2J lN/2\ 

3^ JJ (e - ^f-^^+^ (e + tf-'^+\e' - {z - 1/2)Y-''-^' 

i=l 1=1 

N-2e+l-i 

i {2e + iy {Ae + i) 




(4e + 2i-l) (4e + 2if (4e + 2i + 1) 



(4.4) 



Clearly, this determines a polynomial of maximal degree 2 [A^/2] uniquely. In fact, an 
explicit expression for P{n, N) can immediately be written down using Lagrange inter- 
polation. As it turns out, the resulting expression for P{n, N) is exactly the expression 
covering the last three lines of ( [4.1|) . In view of ([4.2|) , this would finish the proof of the 
Lemma. 

Step 1. D{n,n,N) is a polynomial in n of degree at most N{N + 2). It is standard 
that sums of powers, such as the entries of D{n, n, N), can be expressed using Bernoulli 
numbers. More precisely, we have (cf. [jl^, p. 269ff]), 



n—l m 1 / \ 

s=-n e=o ^ ^ 



where Be denotes the i-th Bernoulli number. Hence, the (i,j)-entry of D{n,n,N) is a 
polynomial in n of degree i + j + 1. Thus, by expanding the determinant D{n,n,N) 
according to the definition of a determinant and determining the degree of each term, 
it follows that D{n,n,N) is a polynomial in n of degree at most X]i<jj<7v(^ + J + 1) = 
2(^+^) + A^ = A^(A^ + 2). 
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From ( [4.5| ) we may read off another property of D{n,n, N), wliicli we record here for 
later use, 

D{-n, -n, N) = {-l)^D{n, n, N). (4.6) 

Step 2. is a factor of D{n, n, N) . From the definition ( p.l|) of the entries of D{n, n, N) 
and ( |4.5|) it is immediate that n divides each entry of D{n,n, N). Hence, divides 
D{n,n,N). 

Step 3. U\=C^ (^^ - ^2)^-2^+1 a factor of D{n, n, N). In view of (|;§, it suffices to 
prove that (n — e)^~2e+i (divides D{n, n, N) for e = 1, 2, . . . , \_N/2\ . In order to do so, we 
claim that for each such e there are N — 2e + l linear combinations of the columns, which 
are themselves linearly independent, that vanish for n = e. Define the coefficients aj by 

e-1 2e-l 

Y[{x - Z) = {X - 6+ l)e-l (X + l)e = ^ ajX^ . (4.7) 
z=—e j=0 

Then we claim that for = 2e, 2e + 1, . . . , we have 

k 

aj_fc+2e-i ■ (column j in D{e, e, N)) = 0. (4.8) 

j=k-2e+l 

To verify (|4.8| ), for i = 1, 2, . . . , we have to prove that 

k e-1 
j=k—2e+l s=—e 

To do so, we interchange sums in (|4.9|), and then use ( |4.7| ), to obtain 

e-1 / 2e-l \ e-1 y e-1 

^ ^.+fc-2e+l J2 a,S^ = n - 

s=—e \ i=0 / s=—e ^ z=—e 

z^a 

for the left-hand side in ( |4.9| ), which clearly vanishes for any integer e > 0. 

Step 4- nlff^(^^ - - 1/2)2)^-2^+1 ^ j^^^^^ of D{n,n,N). In view of (|;§, it 
suffices to prove that (n - (e - l/2))^-2e+i (divides L)(n, n, A^) for e = 1, 2, . . . , [N/2\ . In 
order to do so, we claim, in a manner analogous to Step 3, that for each such e there are 
A^ — 2e + l linear combinations of the columns, which are themselves linearly independent, 
that vanish for = e — 1/2. Define the coefficients aj by 

e-1 2e-l 

Yl (x-2+l/2) = (x-e + 3/2)2e-i= Xl^J^'- (4-10) 

2=-e+l j=0 

Then we claim that for k = 2e,2e + 1, . . . , N we have 

k 

Y aj_fc+2e-i ■ (column j in L)(e- 1/2, e- 1/2, A^)) = 0. (4.11) 

j=k-2e+l 
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To verify ( |4.11| ), for i = 1,2, . . . , N we have to prove that 

E E ^ f ^ ') B.^,-e ((e - l/2Y^' - (-e + 1/2)^+^) = 0, 

j=k-2e+l e=0 ^ ^ 

where we used the right-hand side of formula ( [4 .51) to express the entries of D{n,n,N). 
By a variation of ( [4 .51) , the last equation can be rewritten as 

k e-1 

E «.-A'+2e-i E (^-1/2^^ = 0. (4.12) 

j=fc-2e+l s=-e+l 

Again, to show that ( [4.12| ) is true, we interchange sums, and then use ( [4.10| ), to obtain 

E {{s - xi2r^--^^^ E - i/2)A = E f - \i2r^-^^^^ n (. - .) 

s=-e+l \ jr=0 / s=-e+l \ 2=-e+l 

for the left-hand side in (|4.12| ), which clearly vanishes for any integer e > 0. 

Stt'p 5. Evaluation of the polynomial P{n, N) at n = - \N/2~\ , - \N/2] + 1, . . . , rA^/2] . 
We start by observing that the symmetry relation ( |4.6|) for D{n,n,N) is "inherited" by 
P{n,N). To be precise, we have 

P{-n,N) = P{n,N). (4.13) 
In view of ( [4.13|) it suffices to determine the evaluations of P{n, N) at n = 0, 1, . . . , [A^/2] . 



What we would like to do is, for any e with < e < [A^/2], to set n = e in ( [4.2|) , 
compute D{e,e, N), and then express P{e,N) as the ratio of D{e,e,N) and the right- 
hand side product evaluated at n = e. Unfortunately, this is typically a ratio 0/0 and, 
hence, undetermined. So, we have to first divide both sides of ([4.2|) by the appropriate 
power of {n — e), and only then set n = e. 

This program is easily carried out for e = 0. As we observed in Step 2, each entry 
of D{n,n,N) is divisible by n. Hence, division of both sides of ([4.2|) by n'^ and then 
specializing to 0, transforms (|4.2|) into the equation 



P{0,N) 



deti<jj<7v i'^Bi^j] 



J-J W2J ^2(Af-2i+l) (^j _ ]_y'2)2(Ar-2i+l) ' 

The Hankel determinant of Bernoulli numbers which appears in the right-hand side ex- 



pression can be evaluated by Theorem ^ with n replaced by and a = b = c = d=l. 
Thus, we obtain (^73|). 

The case 1 < e < [A^/2] requires more work. First, we have to "preprocess" the 
determinant D{n,n, N). Define the coefficients aj as before in (|4] 

e-1 2e-l 

Y[{x- z) = {x-e + l)e-i (x + 1), = E 

z=—e j=0 

Then, for k = N, N — 1, . . . ,2e, in this order, add 

fc-i 

a_,_fc+2e-i ■ (column j in D{n, n, N)) 

j=k-2e+l 
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to column k. Thus, by ( |4.8| ), each entry in column k, k = 2e,2e + 1, . . . , N, will be 
divisible by {n — e) after performing these operations. Next, apply the analogous row 
operations. I.e., for k = N, N — 1, . . . ,2e, in this order, add 



fc-i 



aj-k+2e-i ■ (row j in D{n, n, N)) 

j=k-2e+l 



to row k. 

Now we divide (n — e)^~^^+^ on both sides of (f4.2| ), and only then set n = e. As a 
result, from equation ([4.2|) we obtain 



det 



P(e,iV) 



X * 

l<i,j,<N 1 Q Y 



L7V/2J lN/2\ 

""Hie- tf-''+' ll{e + 



-{l- 1/2) 



2\N-2i+l 



i=l 



i=l 



det X det Y 



\_N/2\ 
i=l 



e — I 



[_N/2\ 
1=1 



N-2i+l(^2 



(e^ -{^- 1/2)^) 



2\N-2i+l 



(4.14) 



where X is the (2e — 1) x (2e — 1) matrix X = ( X]s=-e j=i 2e-v ^^"^ where y 
the (iV — 2e + 1) x (A^ — 2e + 1) matrix (yij)ij=i,...,Ar_2e+i whose entries are given by 



IS 



Alternatively, 



Y 



d_ 

dn 



('xi^'^M(^-e + l)e-l(^ + l)e)') 



(4.15) 



The expression on the right-hand side is a certain linear combination of expressions of the 
form 
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In order to compute Hm, we use ( [4.5| ) (and variations thereof) to rewrite it as 



e=o ^ ^ 

i=i ^ ^ 

e-l -1 
s=0 s=—e 

Thus, using the symbohc notation = B^. for the Bernoulh numbers, a combination of 
i ^JBj ] and leads to 

e-l , 

= 2i?^+^- {{B-e + l)e-l (S + l)e)' + i''^' ((^ - e + l)e-l + l)e)') 



s=0 



ds 

s=—e 

2B'+^ {{B-e + l)e-i {B + 1 



as each summand in either sum on the right-hand side vanishes. 
Substituting all this into ( |4.14] ), we arrive at 



p. _ Dj^^ g> 2e - 1) deti<,,,<^_2e+i {2B'+^ {{B - e + {B + 1),)^) 

iN/2\ \_N/2\ ' ^ '> 

fj (e - fj (e + i)''-^'^\e^ - {i - l/lff'^'^' 

1=1 i=l 



which is valid for 1 < e < . The determinant -D(e, e, 2e — 1) is evaluated separately 

in Lemma [T^. The Hankel determinant of (linear combinations of) Bernoulli numbers 
which appears in the right-hand side expression can be evaluated by Theorem ^ with n 
replaced byA^ — 2e + l, a = 6 = e + l and c = d = e. Thus, we obtain 



This completes the proof of the Lemma. 



□ 
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Lemma 13. Let n and m be positive integers. Then the determinant D{n,n — 1, 2m — 1), 
as defined in (|2.1|), is equal to 



-4m+5 



r4m-i .A (2m + 1) (2n-l) 



2 / 2m\ / 2m — 4 



X 



[n-m+i)2„_i (n-ri ^!) ^ ^v^;v^-2 

(2m-l \ /2m-2 
n ^' / I JJ (2ra - z - 2)2i+3 

m (2m - 3) (2m - 1) (n - 2) (n + 1) (n^ - n + 2m + 1) 
(m - 1) (2m + 1) 

6(n - l)n (3),, (I);, (2 - m)h (| + m),. (2 - n)^ (1 + n)h 



+ 



(m + 1) (l);,(|);,(2 + m)4|-m)42 + n);,(3-n)J- ^ 



Likewise, the determinant D{n,n — 1,2m) is equal to 
4 . „x ,,2 /2m -2\ /2m + 2 



(2m + 3) (2n - 1) 



■ 2m \ /2m-2 



X {n- m)m-2 {n + 2)^-2 {n - m - l)2m+2 ^ j |^ (2n - i - 2)2i+3 

(2m - 1) (m + 1) (2m + 1) {n - 2) {n + 1) (n^ - n - 2m) 

m (2m + 3) 

6(n - l)n (3), (f ), (1 - m), (f + m), (2 - n), (1 + n), \ 
(m + 2) ^ (l)4|);,(3 + m)^(|-m)^(2 + n);,(3-n);,y ^ 

Proof. Basically, the proof proceeds in the same way as above. By considerations which 
parallel Steps 1-4 of the previous proof, we deduce that 



D{n,n-1,N) = [n- -] {{n - l)n 



N 

N-1 



2 

r(7V+l)/2l . , . xAr-2i+2x 

n + i^n + t--jin-t)in + t-l)j j-g(n,iV), (4.20) 

where Q{n, N) is a polynomial in n of degree at most 2 \{N + l)/2]. 

Also, Step 5 of the previous proof has a parallel here. Eventually, this yields evaluations 
of Q(n, iV) at n = - \{N + l)/2]+l, - \{N + l)/2] +2, . . . , [(A^ + l)/2l - 1, [(A^ + l)/2l . 
In particular, the symmetry relation which plays the role of (|4.13|) in this new context is 

Q{-n+l,N) =Qin,N), (4.21) 

while the evaluation of determinant D(e, e — 1, 2e — 2) which is needed here (and replaces 
the evaluation of D{e, e, 2e— 1) in this new context; compare ( |4.17] )) is evaluated separately 
in Lemma |1^. 

Unfortunately, this is not good enough. The polynomial Q{n, N) is a polynomial of 
maximal degree 2 [(A^ + l)/2], but by now we have found only 2 [(A^ + l)/2] explicit 
special evaluations of Q{n,N). Hence, we need one more information about Q{n,N). 
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We get this missing piece of information by computing the leading coefficient of Q{n, N). 
This is easily done. By the definition of D{n,n — l,N), given by (p.l|), by the analogue 
of 

n—2 m 1 / \ 

s=-n ^=0 ^ ^ 



and by ( 4.20 ), the leading coefficient is given by 

det if^+jiseven, 
i<i,3<N \ I otherwise. 

It is easy to see, by sorting rows and columns with odd indices to the beginning, that 
this determinant equals the product of the minor consisting of the odd rows and columns 
times the minor consisting of the even rows and columns, explicitly 

det I ^ 1 det 



^<i,3<\N/2\ \2i + 2j -IJ l<i,j<lN/2\ \2i + 2j + 1 

Both of these determinants are easily evaluated by means of the Cauchy determinant 
evaluation (see |^ vol. Ill, p. 311]) 

/ 1 \ ^ Ul<i<j<Ki^^- ^3)iyi-y3) 

l<i,j<K \xi + yj Ul<i,j<Ki^i + 

For the simplifications of the resulting expression for Q{n, N) it turns out to be conve- 
nient to separate the cases = 2m and = 2m — 1. 

This concludes the proof of the Lemma. □ 

Lemma 14. Let n andm be positive integers. Then the determinant D{n, n — 2, 2m + 1) , 
as defined in (|2.1|) , is equal to 

2m+3 



(2m - 2)! (2m + 2)! (n - m - l)2m+i {n-m-2) 



(m - 1)! m! (m + 1)! (m + 4)! (4m + 3)! n(n - 2) 

2m 



-p-r (i + 1)!^ {2n-i- 2)2^+1 / 2 (m + 3) (m + 4) (2m - 1) (2m + 1) x{m, n) 



i=0 



{2m + i + 2)\ y {n-2){n-lfn 
24(2m + 3)(2m + 5)(?2 - m - 3)(n + m + 1) 



(n-3)(n + 1) 

(4)fe (-m)fe (I + m)h (3 - n)h (1 + n)h | 
■^(l);,(5 + m);.(|-m);,(2 + n);,(4-n)^ j' ^ 



where 



x{m, n) = —72m — 204m — 212m — 96m — 16m — 30n — 86mn — 68m n — 16m n 

- 29n2 + Smn^ + 26m'^n'^ + Bm^n^ + 44n^ + 40mn^ + Sm^n^ - lln^ - lOmn^ - 2m^n^. 



22 



M. FULMEK AND C. KRATTENTHALER 



Likewise, the determinant D{n,n — 2,2m) is equal to 

(2m-2)!(2m + 2)!(n-m-l)L+i 

(2m - 3) (m - 1)! ml (m + 1)! (m + 3)! (4m + 1)! n{n - 2) 

{i + 1)!^ {2n-i- 2)2^+1 (m + 2) (m + 3) (2m - 3) (2m - 1) y{m, n) 



i=0 



(2m + 2 + l)! I {n-2){n~lfn 



_ 24m(2m + 3)(n - m - 2)(n + m) ^4 (4)^ (1 - m)^ (| + m)^ (3 - n)h (1 + n)^ 
- 3)(n + 1) ^ (1)^ (4 + m)^ (| - m),, (2 + n);, (4 - n),,^ 

(4.24) 

y(m, n) = 12?ti^ + 44m^ + 48m^ + 16m^ + 16mn + 28m^n + 16m'^n + 4n^ — 20mn^ 

Proof. Basically, the proof proceeds in the same was as in the preceding lemmas. By 
considerations which parallel Steps 1-4 of the proof of Lemma we deduce that 

Af-l / .2\ N-i VN/2\ 



n ((^ - -\) n ((^ - - ) ■ 

(4.25) 

where -R(n, A^) is a polynomial in n of degree at most 2 [(A^)/2'| + 2 if > 2. 



Step 5 of the proof of Lemma has a parallel here, too, which yields evaluations of 
i?(?7,, iV) at n = — [A^/2'l + 1, . . . , 0, 2, . . . , [iV/2'| + 1. In particular, the symmetry relation 
which plays the role of (|4.13| ) in this context is 

R{-n^2,N) = R{n,N). (4.26) 

While the evaluation of the determinant i5(e, e — 2, 2e — 3) which is needed here is done 
separately in Lemma for e > 2, the case e = 2 leads to the 2x2 determinant 
deti<.,,<2 ((-1)'+^' + (-2)^+^) = 4. 

The evaluation of i?(n, A^) at n = 1 requires extra treatment. First, factoring (n — 1) 
out of each column of Diri^n — 2,N) and then setting n = 1 yields the determinant 
deti<jj<Ar (2(5 — 1)*+-^) (using again the symbolic notation B'' = B^). By row and 
column operations, this determinant can be transformed into the form 

2^ det (B'+^-\B-m. 

This determinant agrees with the determinant on the left-hand side of (|5.13|) with n = N , 
a = b = and c = d = 2 (adopting the usual convention (i?+l)_i := 1/B). Unfortunately, 
formula (|5.13| ) does not hold for this choice of parameters. In fact, the evaluation of this 



determinant is rather tedious and therefore given separately in Lemma ^ 

Altogether, we have 2 [iV/2] + 1 evaluations for our polynomial R{n,N) at special 
values of n so far. So we are short of exactly two informations on R{n, N). We can get 
these missing informations by computing the leading coefficient of the polynomial, which 
is done in the same way as in the proof of the preceding lemma, and by exploiting the 
symmetry ( |4.26| ) once again. 
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For the simplifications of tlie resulting expression for R{n, N) it turns out to be again 
convenient to separate the cases N = 2m and N = 2m + 1 . 

This concludes the proof of the Lemma. □ 



2e-l 



Lemma 15. Fore > 1, we have 

D{e, e, 2e-l) = ( £ sA = \{ (4.27) 

\s=—e / i=l 

{The determinant D{a,b,N) was defined in (|2.1|) .) 

Proof. For £ > 1 define the polynomials pe{x) by Pi{x) := {x — \i/2'] + 1)^. Furthermore, 

define coefficients a^^k, = 1, 2, by pe{x) = J2k=i '^i,kX^ ■ Now, for £ = 2e — 1, 2e — 

2, . . . , 2, in this order, we replace column q of the matrix underlying -D(e, e, 2e — 1) by the 
linear combination of columns 'Yl,k=i (^i,kCk, and afterwards do the same sort of replacement 
in the rows. These operations yield a matrix with (z, j)-entry equal to 

e-1 e-1 

^ p.{s)p,{s) =Y.{s- rV2l + 1). {s - \j/2] + 1),. (4.28) 

s=—e s=—e 

The next sequence of operations will turn this matrix into upper triangular form, so 
that the determinant is easily obtained by forming the product of the diagonal entries. 

To begin with, it should be noted that for j = 2e — l the sum ([4.28|) consists of just the 
term corresponding to s = — e. We subtract pj{—e) /p2e-i{—e) times column 2e — 1 from 
column j, j = 1, 2, . . . , 2e — 2. The previous observation tells us that these operations 
have the effect that the entries that are in one of the columns j = 1, 2, . . . , 2e — 2 become 

e-1 e-1 

E p^i')pM = E - rv2i + 1). {s - \jm + 1)„ (4.29) 

s=— e+l s=— e+1 

i.e., the summand corresponding to s = — e has been eliminated. In particular, all the 
entries in the last row, except the rightmost entry, of course, are 0. 

Next consider column 2e — 2. After the above column operations, the sum ( 4.29| ) with 



j = 2e — 2 which defines the entries collapses to just the term corresponding to s = e — 1. 
We subtract pj(e — l)/p2e-2(e — 1) times column 2e — 2 from column j, j = 1, 2, . . . , 2e — 3. 
The previous observation tells us that these operations have the effect that now the entries 
that are in one of the columns j = 1, 2, . . . , 2e — 3 become 

e-2 e-2 

E Ms)p,{s) = E - + 1)^ - r^V2l + 1)„ (4.30) 

s=— e+l s=—e+l 

i.e., the summand corresponding to s = e — 1 has been eliminated as well. In particular, 
all the entries in the next-to-last row, except the two rightmost entries, of course, are 0. 

If we continue in the same manner, then eventually we arrive at an upper triangular 
matrix whose i-th diagonal entry is il"^. Thus, the result ( ^.27| ) follows. □ 

Lemma 16. For e > 2 we have 

D{e, e - 1, 2e - 2) = ^^^det^_^ | > ] s^^^ ] = \ \ (4.31) 
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{The determinant D{a,b,N) was defined in ( |2.1| ).) 

Proof. This identity can be established in essentially the same manner as the previous 
lemma. The basic difference is that one has to replace the polynomials pe{x) in the 
previous proof by the polynomials 




if £= 1, 
-[i/2\+l)e ife>l. 



Everything else is completely analogous. We leave the details to the reader. □ 
Lemma 17. For e > 3 we have 

Die,e-2,2e-3) = ^^det(Y.'^'"') = U^^-'- (4-32) 

\s=—e / i=l 

{The determinant D{a,b,N) was defined in (|2.1|) .) 

Proof. This identity can be established in essentially the same manner as the previous 
lemmas. The basic difference is that one has to replace the polynomials in the previous 
proofs by the polynomials 

{X if £ = 1, 

x{x + 1) if £ = 2, 

{x - \i/2] + 2)e if £ > 2. 

Everything else is completely analogous. We leave the details to the reader. □ 



5. Orthogonal polynomials, continued fractions, and Hankel 
determinants of bernoulli numbers 

In this section we review some facts about the interrelations between orthogonal poly- 
nomials, continued fractions, and Hankel determinants. Good sources for information 



about these topics are |T8|, ^ |7|' 0' H' 
To begin with, we recall Favard's Theorem. 

Theorem 18. (Cf. ||48| , Theoreme 9 on p. 1-4] or ||4^, Theorem 50.1]). Let {pn{x))n>o be 
a sequence ofmonic polynomials, the polynomial Pn{x) having degree n. Then the sequence 
{Pn{x)) is (formally) orthogonal with respect to a linear functional L , i.e., L{pn{x)pm{x)) = 
SmnCn for somc scqucncc (c„)„>o of nonzero numbers, with 6m,n denoting the Kronecker 
delta (i.e., 6m,n = 1 if m = n and 6m,n = otherwise) if and only if there exist sequences 
{o-n)n>i o,nd {bn)n>i, with bn ^ for all n > 1, such that the three-term recurrence 



Pn+l{x) = {an + X)pn{x) - bnP„~l{x) , for U > 1, 

holds, with initial conditions po{x) = 1 and pi{x) = x + Oq. 



(5.1) 



It is a simple fact that, given a linear functional, the corresponding orthogonal polyno- 
mials can be expressed in form of certain determinants. 
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Theorem 19. (Cf. [|9|, Theorem 7.15, (7.2.41)] or ^ (36)]). Let L be a linear functional 



defined on polynomials with moments /i„ 
of orthogonal polynomials is given by 



L{x^). Then the corresponding set {pn{x))n>o 









fJ'2 














/in 




Pn{x) = (i„Met 


fJ'2 




















f^2n-l 




\ 1 


X 






x" / 



(5.2) 



where dn = det 



0<i,j<n-l{fJ'i+j ) ■ 

The next theorem addresses the relation between orthogonal polynomials and continued 
fractions, the link being the generating function of the moments. 

Theorem 20. (Cf. p]. Theorem 51.1] or |4|, Proposition 1, (7), on p. V-5]). Let 

{Pn{x))n>o be a sequence of monic polynomials, the polynomial Pn{x) having degree n, 
which is orthogonal with respect to some functional L. Let 

Vn+l{x) = (a„ + X)pnix) - bnPn-lix) (5.3) 

be the corresponding three-term recurrence which is guaranteed by Favard's theorem. Then 
the generating function for the moments /i„ = L(x") satisfies 



n=0 



X" 



1 + QqX 



bix 



(5.4) 



1 + aix 



b2X' 



1 + a2X 



We remark that a continued fraction of the type as in ( |5.4| ) is called a J-fraction. 
The next theorem addresses the relation between J-fractions and Hankel determinants. 



Theorem 21. (Cf. H, Theorem 51.1] or ||, Corollaire 6, (19), on p. IV-17]). Let 

{f-n)n>o be a sequence of numbers whose generating function, when written in terms of 
a J-fraction, is given by (|5.4|). Then the Hankel determinant deto<ij<n-i(/ii+j) equals 

nun-lhn-2 i2 I. 

The family of orthogonal polynomials that is of interest to us is the continuous Hahn 
polynomials, which were first studied by Atakishiyev and Suslov |^ and Askey p. The 
theorem below lists the relevant facts about these polynomials. 

Theorem 22. (Cf. Section 1.4]) . The continuous Hahn polynomials {pnid, b, c, d; x))„>o 
are the monic polynomials defined as the terminating hypergeometric series 



Pn{a,b, c,d;x) 



^)" (a + c)n (a + d)n 



-n,n-\-a-\-b-\-c-\-d — l,a-\-xy 
a -\- c,a -\- d 



;1 



{a + b + c + d + n — 1)„ 
These polynomials satisfy the three-term recurrence 

Pn+i{a, b, c, d; x) = {x - A„(a, 6, c, ci))p„(a, b, c, d; x) - Bn{a, b, c, d)pn-i{a, b, c, d; x) 



(5.5) 



(5.6) 
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where 

A ( h A\- ^ n (6 + c + n-l) ih + d + n-l) 

/i„^a,o,c,aj - V i (^«+ (^^^^^^^^ 2n- 2) (a + 6 + c + d + 2n-l) + 

[1 — a — h — c — d — n) (a + c + ra) {a + d + n) 
{-l + a + h + c + d + 2n) {a + h + c + d + 2n) 

n (a + c + n — 1) (h + c + n — 1) {a + d + n — 1) 



(5.7) 



and 
Bn{a, b, c, d) 



(a + 6 + c + o? + 2n-3) [a + h + c + d + 2n - 2) 



2 



(6 + rf + n-l)(a + 6 + c + rf + n-2) 
^ ia + h + c + d + 2n-l) ' ^ ' 



They are orthogonal with respect to the functional L which is given by 

L{p{x)) = — T{a + xy/^) r{b + x^/^) r(c - xV^) r{d - x^/^) p{x) dx . 

(5.9) 

In this integral, if necessary, the path of integration has to be deformed, so that it separates 
the increasing sequences ((a+A;)-\/— l)^,^^ and ((6+A;)-\/— l)^.^^ of poles from the decreasing 

sequences ( — (c + A;)-\/— l)^>o + k)\/—l)^^^ of poles (see 0/ 

Explicitly, the orthogonality relation is 



L{pm{a, b, c, d] x)pn{a, b, c, d; x)) 

r(n + a + c)T{n + a + d)T{n + b + c)T{n + b + d) n\ 



{a + b + c + d + n — 1)„ (a + fe + c + d + n — l)„+ir(n + a + b + c + d— 1) 

(5.10) 

In particular, 

(g + c-l)! (fe + c-1)! (g + ^-l)! {b + d-l)\ 
^^^^ - {a + b + c + d-l)\ • ^^-^^^ 

Remark. (1) The reader should be aware that Theorem ^ is formulated for the monic 
form of the continuous Hahn polynomials, so that our polynomials are not the polynomials 
which are denoted by p„(x; a, 6, c, d) in Section 1.4] but those which are denoted there, 
slightly confusingly, by Pn{x). 

(2) Although the definition ( p. 5] ) does not show it, the continuous Hahn polynomials are 
symmetric in a and 6, and in c and d, because the orthogonality measure (|5.9|) has these 
symmetries. In addition, there is another symmetry when the roles of a, b are interchanged 
with the roles of c, d. Namely, we have 

p„(a, 6, c, d] x) = (— l)"'j9„(c, d, a, b; —x). (5-12) 

This follows from the fact that, when performing the substitution x —x in the integral 
( p.9| ), we obtain almost the same integrand, the only differences being that the roles of 
a, b and c, d are interchanged, and that p{x) is replaced by p{—x). The path of integration 
also does not change, at least as long as all of a, b, c, d are positive. This is, however, 
sufficient to conclude that ( p.l2| ) holds for all a, b, c, d because, for fixed n, both sides of 
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( p.l2| ) are rational of bounded degree in a, 6, c, d. (The same conclusion could also be less 
elegantly derived by applying some 3F2-transformation formulas.) 



Now, by combining Theorems [2^, and |22| we are able to derive without difficulty 
the determinant evaluation that we need in Step 5 of the proof of Lemma O and the 



analogous places in Lemmas |T3| and |T4| (which, in turn, are essential for the proofs of 
Theorems |1]-§D. 

Theorem 23. For positive integers a, b and nonnegative integers c, d there holds 



det (5*+^' {B + {B + i-B + i-B + l)d„i) 

l<i,j<n 

^_i)n(n~i)/2 + (b + c-iy. {a + d-i)\ {b + d-iy.y 



n-1 



X 



n 



(a + b + c + d-l)\ 



i (a + c + i-l) (b + c + i-l) (a + d + i - I] 



\\{a + b + c + d + 2i-3) {a + b + c + d + 2i - 2) 

(b + d + i-l) {a + b + c + d + i-2) 
(a + b + c + d + 2i-l) 



, (5.13) 



where we have again used the symbolic notation B^ = B^. (In case that c or d should be 
zero, we have to interpret {—B + l)_i as l/{—B).) 

Proof. We claim that the moments for the continuous Hahn polynomials, i.e., for the 
linear functional as defined in (15.91), are 



/i„ = (B/V-irB' {B + l)a-l {B + l),_i i-B + l)e-l i-B + 



(5.14) 



Let us suppose for the moment that this is already established. Then, by Theorem 20 



and (|5.6| ) of Theorem |2^, we have 



[B/^/^rB^ {B + l)a-i {B + 1)5-1 i-B + l)e-i i-B + l),_ia;" 



n>0 



52 {B + l)„_i {B + 1)5-1 i-B + l),_i i-B + l)d_i 



1 + Ao(a, b, c, d)x 



Bi{a, b, c, d)x 



1 + Ai{a, b, c, d)x 



B2{a, b, c, d)x 



1 + ^2(0, 6, c, d)x 



But, by Theorem this immediately implies the truth of ( [5. 131) 
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It remains to verify ( [5.14| ). Using ( |5.9| ), this is a rather straightforward computation: 
1 



— I r(a + XA 



-i)r(6 + xA 



-1 r c 



1 

2^ 



CO 

ooV— 1 
— oov'— 1 

1 



r(a + z) T{b + z) r(c - z) V{d - z) 



-l)T{d- xV^) dx 



dz 



2^T{^ 



-1 



n+1 



{z)a {z), i-z), i-z), {T{z) T{-z)f z^ dz 



2tt( 



— oov — 1 

OO-v/— 1 



-Z c [-Z d 



2tt( 



■I) 



n+1 



oov — 1 

OO-v/— 1 



TT 



z sm vrz 



{z + (2; + (-2; + (-Z + 1) 



d-1 



CXDV — 1 



vr 



sm TTZ 



(5.15) 



In the third line we used the relation T{a + z) = {z)a T{z) for a G N (see, e.g., |T^, 1.2(2)]), 
and in the fourth line we used the formula T{z) T{—z) = — ^^.^^^ (see [12, 1.2(5)]). The 



reader should notice that, because of the convention regarding the path of integration 
in ( ^.9|) , in case that c or c? are zero the path of integration in the last line of ( ^.15|) is 
deformed so that it crosses the real axis between the singularities z = —1 and z = 0. 

To finish the calculation, we appeal to the following integral representation of Bernoulli 
numbers (see [p3| , p. 75], and let a — > O"*" in the relevant identity in the middle of the 
page) 

fOO-v/— 1 



1 



27r. 



4 



00\/—l 



TT 



smiTZ 



dz. 



(5.16) 



(If = or z/ = 1 then the path of integration is indented so that it avoids the singularity 
z = 0, passing it on the negative side, but to the right of the singularity z = —1.) If we 
use this formula in (|5.15|) , we obtain ( |5.14D immediately. □ 

At least two special cases of Theorem |2^ have explicitly appeared in the literature 
before. The case a = b=l, c = d = appears for example in |]l|, (3.1)]. The case 
a = b = c = d = l appears in ||30| , App. A. 5, p. 322]. 

In the proof of Theorems | and |^ (to be precise, in the proof of Lemma ^5]), we make 
use of a rather recent result on (formal) orthogonal polynomials, due to Leclerc [^, 
Theorem 1]. 

Theorem 24. For an arbitrary sequence of numbers (/i„)„>o let {Pn{x))n>o be the se- 
quence of polynomials defined by 

/in \ 

A^n fJ'n+1 
/in+1 /in+2 





/ /io 


Ail 


/i2 










Pn{x) := det 






A^n 




l^n-l 




A^n+l 



\ 1 



X 



X 



n-1 



A'2n-1 

X^' 



(5.17) 
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(These are, up to normalization, the orthogonal polynomials with respect to the linear 
functional with moments (/in); compare Theorem^^.) Furthermore, let {Qnix))n>o be the 
sequence of polynomials defined by 



k=o ^ ^ 



—X 



\n—k 



Then, for arbitrary integers I > 1, m > 1, there holds 



( Pii 



det 



x] 



Pl+l{x) ... Pl+m-lix)\ 



P/(x) P'{x) ... Pi 



\Pt''\x) Pl--'\x) 



+m— 1 



X] 



Ci.m. det 



/ Qm{x) Qm+l{x) 
Qm+l{x) Qm+2{x) 

\Qm+l-l{x) Qm+l{x) 



where 



m— 1 



Lm 



-1)'" Yl A;! det 



k=l 



\P-k+l-l f^k+l 



Qm+l-l{x) \ 
Qm+l{x) 



Q 



m+2l-2 



f-k+l-l \ 
f^k+l 

f^2k+2l-2/ 



(5.18) 



(5.19) 



6. Auxiliary lemmas 



Lemma 25. For any integer n > 1, we have 
det (B'+^-\B-l?) 

n 

= (_i)"("-i)/2g-n(i ^2n + 5n^ + An^ + n") JJ 



i=l 



z(z + l)4(i + 2) 



(2i + l)(2i + 2)2(2i + 3) 



(6.1) 



(again, using the symbolic notation = Bj.). 

Proof. The determinant in (6J.) is equal to the determinant in ( 5.13|) with a = b = and 
c = d = 2 (again, with the convention that (i?+l)-i is interpreted as 1/B). Unfortunately, 
for this choice of parameters, formula (|5.13|) is not valid. However, the determinant in 
( p.l|) is very close to the determinant in ( ^.13| ) with a = b = 2 and c = d = 0. In fact, 
because of the well-known property of Bernoulli numbers that B2k+i = for all positive 
integers k, multiplication of all even numbered rows and columns of the determinant in 
( p.l3| ) by —1 (which does not change the value of the determinant) with this choice of 
parameters transforms the latter into the determinant := deti<jj<„(Aj+j„2); where 

Xo = B\B + iy = B\B -lY -2, 
Ai = -B\B + 1)2 = B\B - 1)2 + 1, 

A„ = (-1)"5"(5 + 1)2 = B'\B - 1)2 for n>2. (6.2) 
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So, because of the deviating definitions of Aq and Ai, the only difference between A^t 
and the determinant in ( |6.1| ) is in the top-left entry and its right and bottom neighbour. 
Linearity of the determinant A„ in the first row and column then implies that 



det {B'+^-\B - If) = A„ - Aii'2}^ii.2} + 2A^i}^^i} + 2/\\l^'^^\ 

l<i,j<n 



(6.3) 



Here, ^{n,i2,...};{jij2...} denotes the minor of A with rows ii, • • • and columns j'l, j2 • • • 
deleted. (Empty minors are defined to be zero). 

Now observe that from Theorem with a = b = 2 and c = d = Owe obtain immediately 
that 



A. 



Likewise, from Theorem 23 with a = b 



i=l 



(22)!(2z + l)! ■ 



(6.4) 



1 and c = d = 2 we have 

n-l 



1=1 



(,_!)! (, + 1)14 (, + 3)1 
(2i + 2)!(2z + 3)! 



(6.5) 



Our next observation is that aIi^'^'^^ is, essentially, the coefficient of x in the continuous 



Hahn polynomial pn-i{0, 1, 2, 2; x). To make a more precise statement, consider 
/ifc = Afc+i, k = 0,1, . . . . Then, obviously, A^^'^'^^'^ equals 

(—1)" det (Ai+,_i) ■ (coefficient of x in p„_i(x)). 

l<jj'<n— 1 



with 



On the other hand. Theorem |I9] says that, with this choice of the /i^, the polynomials 
Pn{x) are orthogonal with associated moments A^+i, k = 0, 1, ... . By comparing 
and ( |5.14 ) with a = b = 2, c = 0, and d = 1, we see that we must have 

p„(x) = (- v^)>„(2, 2, 0, 1; -x/V^). 

Therefore, if we remember (|5.12| ), we obtain that Ai^^'^^^ equals 



'-in 



-1 



>n-l 



det (Ai+,_i) • (coefficient of x in p„_i(0, 1, 2, 2; x/a/^)) . (6.6) 

l<ij<n— 1 



The Hankel determinant in this expression can be evaluated by using Theorem ^ with 
a = b = 2, c = 0, d = 1, and n replaced hj n — 1. By substituting the result in (|6.6|) and 
by using the definition ( |5.5| ) of continuous Hahn polynomials with a = 0,b=l,c = d = 2, 
we obtain 



A^^^-'f^^ = (_i)"("+i)/2 



nl 



{n + 3)„,_i 



n-l 

n 



(22 + l)!(2i + 2)! 



(, _ 1)1 ,12 (, + 1)12 (, + 2)!\ ^ (1 - n)u (n + 3)^ 



k=l 



k{k + iy? 



(6.7) 



The remaining minor 

^{1,2};{1,2} 

requires additional work. We employ Theorem ^ with 
I = n — 2, m = 2, X = 0, and /i^ = Afc+2, k = 0,1, . . . . With this choice of parameters, the 
determinant on the right-hand side of ( p.l8| ) is precisely our remaining minor An ' ''^ ' \ 
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The single determinant occurring in the definition ( |5.19| ) of Cn-2,2 can be evaluated by 



using Theorem p3| with a = 6 = l,c = (i = 2, and n replaced by n — 2, so that we have 



n-1 



+ (2 + 3)! 
(2z + 2)!(2i + 3)! 



(6.8) 



The determinant on the left-hand side ( |5.18|) is a 2 x 2-determinant whose entries are 
the constant term and the coefficient of x, respectively, of Pn-2{,x) and Pn-i{x). The 
polynomials Pn(x), n = 0,1,..., on the other hand, are orthogonal polynomials with 
associated moments Xk+2, k = 0,1, . . . (compare ( p.l?| ) and Theorem |19D. By comparing 
the definition ( |6.2| ) of the Aj's with ( |5.14|) , it is seen that the polynomials Pn{x) must agree, 
up to normalization, with the continuous Hahn polynomials with parameters a = b = 1 
and c = d = 2. To be precise, we have 

Pn{x) = {V^r det (A,+,)p„(l,l,2,2;a;/v^). 

l<t,j<n 

Clearly, the Hankel determinant on the right-hand side can be once again evaluated by 
means of Theorem |2^ with a = b = 1 and c = d = 2. In summary, from ( ^.18|) with the 
above choice of parameters we infer 



A 



{1,2};{1,2} ^ ^ ^■^n{n-l)/2 (3)n-2 (3)n-l 



^n-2 

n 



{n + 3)n-2 {n + 4)„_i (2^ + 2)! {2t + 3)! 

X (Co,n-2Cl,n-l ~ Co,n-lCl,n-2) ; (6-9) 



where co,n and ci,„ denote the coefficient of x° and x^, respectively, in 

{-n)k {n + 5)fc (1 + x)k 



E 

k=0 



k^- (3)^ 



Using hypergeometric notation ( p.6|) , the first of these two, Co,n, can be expressed as 



[1 + n) (2 + n) (3 + n) (4 + n) 



-Fi 



—2 — n,3 + n 



;1 



-l + (n + 2)(n + 3) 



The 2-Fi-series can be evaluated by means of the Chu-Vandermonde summation formula 
(see H, (1.7.7); Appendix (III.4)]) 



2^1 



a, -N_ 

1 -L 



(c - a)N 



(c) 



N 



(6.10) 



where iV is a nonnegative integer. This yields 



Co,', 



,,w , if n is even, 

n+l)(n+4) ' 



(n+2)(n+3) 



if n odd. 



(6.11) 
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By combining ( |6.7| ), ( |6.9| ) and ( |6.11| ), we get, after a considerable amount of simplification, 

^ ^ J-J- (2z)!(22 + l)! 



'n-2 



-1)^ (n + + 3)! (1 + (-l)"(n + 1)) 

n-2 



{n - 1)! ^ (A; + 2)!(A; + 3)!(n- A;-2)! j + 



1 \ 



By Lemma p6|, with n replaced by n — 1, this expression reduces to 



2A{i};{2} _ A{i,2};{i,2} ^ + 1) + 2) 



(n + 2)! A (z-l)!^!" (z + 1)! 
{n-iy.l\ (20!(2z + l)! 



(6.12) 



Substituting ( |6.4|) , ( |6.5|) , and (|6.12| ) in ( |6.3|) , and simplifying the resulting expression, we 
eventually arrive at ( |6.1|) . □ 

Lemma 26. For n > 0, we have 



(n + 3)! 



n-l 
fc=0 



(-1)^ {n + k + 4)! (1 - + 2)) 



n— 1 
fc=0 



(A; + l)(A; + 2)!2(n-fc-l)! 



(A; + 2)!(fc + 3)!(n- A; - 1)! + 



k 

= (-!)> + 2) -2. (6.13) 



Proof. We shall treat the two sums (in the first and the second line in (|6.13|) , respectively) 
separately. 

First we consider the sum in the first line of ( 3.13| ). We replace the term l/(/c + 1) by 
x'^ dx, interchange summation and integration, and write the sum in hypergeometric 
notation ( |2.6| ). This gives 

1 - (-l)"(n + 2) fn + 2 



(n + l)s 



x^ 



3-r 2 



n + 3, —n — 1 



X 



X^ 



X 



Using the transformation formula (see (1.7.1.3), sum reversed at the right-hand side]) 



2^1 



a, -N 



jv (c-a)jv 
(c)iv 



-N, l-c-N_ z - 1 

1 + a- c- N' z 



(6.14) 



(where iV is a nonnegative integer) with a = n + 3, A^ = n+1, c=l, and z = x, this is 
transformed into 



1 - (-l)"(n + 2) 



„1 / n+l 

-^0 V fc=o 



A;! (A; + l)! 



(6.15) 



x^ 



X 
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Recall that for 3?(q;) > —1 and 3?(/3) > —1 we have the following integral representation 
for the Euler beta function, 



^ ^ T{a + p + 2) 



(6.16) 



Use of this identity in ( |6.15| ) wherever possible (i.e., it is applied to the summands with 
k = 0, . . . ,n — 1) yields after some simplification the expression 



1 - (-l)"(n + 2) 
in + 1)3 



k=0 



i-n-2)l 



^ ^•Y (n + l)2((n + 3)-(n + 2)(l-xr) ^ (n + 2) {{1 - xT^' - l) ' ^ ^ ^^^^^^ 



X 



x^ 



We would like to write the sum in the first line as a hypergeometric series. Unfortunately, 
this cannot be done by just straighforwardly extending the summation over all nonnegative 
k because of the term {—n)k in the denominator, which is for k = n + 1. The way to 
overcome this problem is to rewrite the sum as a limit. 



fc=0 



-n 



2)1 



[-n)k k\ 



lim 



—n — 2,n — 2 



-n — e 



;1 



-n-2f 



n+2 



(n + 1)\ {-n - e)n+i (n + 2)\ {-n - e)n+2 



The 2-^1-series can be evaluated by means of the Chu-Vandermonde summation formula 
( p.lU| ). Substituing the result in ( |6.17D , we obtain 



1 - (-l)"(n + 2)(^ / _ / (2 - eU, _ (^ + 2) (n + 2)! (1 - e) - (n + 2)! 



n + 1 



-n - e)n+2 



' /^( n + l)2((n + 3) - (n + 2) (1 -x)") ^ (n + 2) ((1 - - 1) ' ^ 



x 



x" 



Of course, the limit can be computed by de I'Hospital's rule, so that we get 



1 - (-l)"(n + 2) 



;^((- + 3)!E^-(- + 2)!(n + 2)j 



(n + 1)3 \ ' ■ n\ 

'1 /(n + l)2((n + 3) - (n + 2)(l-x)'^) (n + 2) ((1 - x)"+i - 1) ' 



+ 



X 



+ 



dx 



x^ 
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Now, let us turn to the integral. Expanding the integrand by the binomial theorem and 
simplifying the result, we obtain 



1 - (-l)'^(n + 2) 
(^ + 1)3 



-1) 



n+l 



[n 



+ 



-l)"(n + 2) 



(n + 1); 



[n 



„1 /n-l 

+ 2)/ E 



n+2 ^ \ 

i)(^ + 2f + (^ + i)3Er-T 

fc=l / 



'-"-^>^^'-'*^ + f'-"-^>'-^^.^ I d.. (6.18) 

yk -\- 2)1 



Now the integration can be performed without any difficulty. We rewrite the resulting 
sum over k in the second line limit, 



E 

k=0 



lim 

£^0 



-n - l)k+2 - {k + 2){-n - 2)k+3 



(k + l) {k + 2)\ 



n+l n+l 



—n, e 



;1 



(n + l)(n + 2) (n + l)(n + 2) ^ 

H 2-ri 



-n, e 



Again, Chu-Vandermonde summation ( |6.1(]| ) can be applied to evaluate the two 2-^1-series. 
If the result is substituted in ( |6.18|) , we get 



1 - (-l)"(n + 2) 



in + l] 



- in 



n+2 ^ \ 

+ 1) {n + 2f + {n + l),Y,-— \ 

k=i J 



+ 



1 - (-l)"(n + 2) f{n + l)\-{2-e 



n + 3 



lim 

£^0 



(n + l)\e 



+ {n + 2) 



n\-(l- 



nl e 



Using de I'Hospital's rule once more, we can compute the limit and obtain after some 
simplification the expression 



:i-M)>+2))^ 



2 (-1)" + 2^3 + 11^2 + 19n + 11 



,k=0 



(n + l)(n + 2)(n + 3) 



(6.19) 



Now we turn our attention to the double sum in the second line of (|6.13|) . Analogously 
to before, we replace the term l/(j + 1) by dx. This enables us to evaluate the inner 

harmonic sum Yl^j=Q ^/U + 1) Iq ~ ^''~^^) / ~ ^) We substitute this in the double 
sum in the second line of ( |6.13| ). Using hypergeometric notation, the result is 



-1)^^(71 + 2) 



1 — X 



n + 3, —n — 1 



— 2F1 

X 



n + 3, —n — 1 



X 



H — I dx. 

X . 



The first hypergeometric series can simply be computed by Chu-Vandermonde summation 
( p.lOp . To the second hypergeometric series we apply the transformation formula 
(1.8.10), terminating form]) 



see 



46 



iFr 



-N 



a)N 



Fi 



a, 
1 + a 



-N 



N' 
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where iV is a nonnegative integer. These operations yield 

2F1 







'1 



X) 



X[\ — x) 



n + 3, —n — 1 



;1 



X 



+ 



n + 2 



X 



dx. 



(6.20) 



Now we would like to apply the Euler beta integral formula ( |6.16| ) once more. However, 
this is not possible just straightforwardly, because the beta integral on the left-hand side 
of (|6.16 ) is not defined for (3 = —1. In order to overcome this problem, we first rewrite 
the term l/x(l — x) (which appears in the second term of the integrand in ( |6.20| )) as 
1/x + 1/(1 — x), and then replace all occurrences of 1/x by \mie^Q+ 1 so that ( |6.2CI| ) 
becomes 



'n+l 



n+l 



X 



+ Y: ^-^^±34z!^x^ (1 - xr^ + i-mn + 2)x-^ 



fc=l 



dx. 



Next we interchange limit and integration, and apply the Euler beta integral formula 
( p.l6| ) wherever possible. In the result, the first sum can easily be evaluated by Chu- 
Vandermonde summation (|6.10|) . Subsequently, we compute the limit by using de I'Ho- 
spital's rule again. This yields the expression 



(n _ 71+1 , \ n+l 

k=0 k=l / k=l 



-n 



Replacing once more the term 1/fc in the rightmost sum by Jq x^ ^ dx, we obtain the 



expression 

f-1 



n 



+ 2) (EiiT+Erii)+/'; 

\fc=0 k=l / -^^ 











^ -I ^ 


- 2F1 







n + 3, —n — 1 



X 



dx. 



To the 2-^1-series we apply the transformation formula ( |S.14| ). In the result, we replace 
any occurrence of x"^ by lim£_»o+ x™"'"'^, so that we arrive at the expression 



•-ir(n 



+2) (Ei^+Ei^) 

\fc=0 k=l / 



+ 



n+l 



lim 



-x'-^ + (-1' 



n+l I 



n 



+ 2)E 



k=0 



{-n-l)l 
k\ {k + l)\ 



{x - Ifx 



k ^n—k+e 



dx. 



Again, we interchange limit and integration and apply the Euler beta integral formula 
( p.l(j| ) once more. Writing the result in hypergeometric notation, we obtain 



(n n+l 
Ei^-Ei^ 
A;=0 k=l 



+ lim ( -- + (-1)"+!^^ 



1- 2F1 



-n — 2, —n — 2 
—n — 1 — £ 



;1 
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The hypergeometric series can be computed by Chu-Vandermonde summation ( |6.10| ). 
Once more we make use of de I'Hospital's rule for the limit and obtain after some simpli- 
fications the expression 



l + (-l)"(n2 + 3n + 3) 



l-i-inn + 2))J2T—T- (6.21) 



n + 2 ^ \ y \ k + 1 



k=0 



When adding together ( |6.19| ), the result of our computation for the sum in the first line 
of ( |6.13| ), and ( |6.21|) , the result of our computation for the double sum in the second line 



of ( |6.13| ), the harmonic sums cancel, and it is easy to verify that, magically, the remaining 
terms simphfy to the right-hand side of ( |6.13| ). □ 



Remark. Peter Paule demonstrated to us, that the identity ( |6.13| ) can also be proved 



algorithmically. Clearly, the Gosper-Zeilberger algorithm ^ finds a recurrence 

for the sum in the first line of ( |6.13| ). Carsten Schneider's extension of Karr's algorithm 
[19|, implemented by Schneider, finds a recurrence for the double sum in the second line 
of ( |6.13| ). Finally, Mallinger's Mathematica package GeneratingFunctions or Salvy 
and Zimmermann's Maple package gf un can be used to combine these two recurrences 



into one, a recurrence of order 10. It is then routine to check (preferably on the computer) 
that the right-hand side of ( |6.13 ) satisfies this same recurrence. However, in the present 



implementation, these algorithms are not able to find the explicit evaluations, in terms of 
harmonic numbers, of the sums in the first and second line of (|6.13|) , given in (|6.19| ) and 
( p.21| ), respectively. 



7. Evaluations of Hankel determinants featuring Bernoulli polynomials 

There are several theorems hidden in the body of this paper. Among these are evalua- 
tions of Hankel determinants of Bernoulli polynomials evaluated at special values. Recall 
that the Z-th Bernoulli polynomial is defined by 

I 



k=0 ^ ^ 



The Hankel determinants of Bernoulli polynomials of which we are talking are special 
cases of the determinant B{N\ x) given by 

B{N-x) := det (5,+,(x)). (7.1) 

\<%,J<N 

This is in a fundamental way different from the Hankel determinant 

which has been considered earlier (see Sec. 5]). (Note that the difference is that, in 



the latter determinant, indices start already with 0.) As is not difficult to see (cf. 



p. 419] or Lemma 15]), the latter determinant does in fact not depend on x (i.e., the 
powers of x cancel in the expansion of the determinant), so that its value is equal to its 
value at X = 0, which, in turn, is given by Theorem ^ with a = 6 = l, c = (i = 0. This is 
in sharp contrast to the Hankel determinant ( [7.1| ), where the powers of x do not cancel, 
so that ([7.1|) is a nontrivial polynomial in x. As such, the evaluation of the determinant 



( [7.1D is much more difficult. Below, we provide evaluations of ( [7.1D for x = —1, x = — 1/2 
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and X = 1/2. Needless to say that the evaluation in the special case x = (and as well 
in the special case x = 1) is given by Theorem with a = b = c = d=l. 

First of all, in the proof of Lemma 0, we observed that (in symbolic notation = 3^) 

det iiB - ly+'J) = det {B'^'~\B - if). 

^<i,j<N ^ ' l<i,j<N ^ ' 



The determinant on the right-hand side was then evaluated in Lemma |2^. The linear 
combination [B — 1)' of Bernoulli numbers is nothing else but i?;(— 1), the /-th Bernoulli 
polynomial evaluated at —1. Thus, we obtain the following corollary. 

Corollary 27. Let N be a positive integer. Then there holds 
BiN; -1) = ^ det ^ (5,+, + (-ly^^t + j)) 



N 



.lf(N-i)/2Q-N^-^ + 2iV + 5iV2 + 4iV3 + N') Y[ 



i=l 



i{i + l)\i + 2) 



(2z + l)(2i + 2)2(2i + 3) 



N-i 



(7.2) 



Our next evaluation results from the determinant evaluation in Lemma [13. 
Theorem 28. Let m be a positive integer. Then for 



B{N- -1/2) 



det ((2^-'-^ 



1)5,+, - (-1/2) 



i+j-li 



+ J)) 



there hold 
B{2m- l;-l/2) 



\m— 1 



X 



3 + 8m + 



{2m -ir (^■^^'^!)'(^^=l(2^-l)!)^ 
26(m-i) (m-i)!6(n;rr'^!) 

(m - 1) (2m + 1) 



3m (m + 1) (2m — 3) (2m 



3, |,2 - m, I + m_ ^ 
f,2 + m, f - m ' 



(7.3) 



and 

B{2m- -1/2) = 
X I 1 + 8m — 



8m (2m + 3) 
3 (m + 1) (m + 2) (2m - 1) (2m + 1) 



L 2 



3 

, 3 + m, I — m 



3, 1, 1 — m, + m 



;1 



(7.4) 



Proof. Consider the determinant -D(n, n— 1, N) (see ( |2.1| ) for definition). Factor 2(n — 1/2) 
out of each column of Z)(n, n — 1, N), and then set n = 1/2. By the appropriate variant 
of ( [4. 51 ) and de I'Hospital's rule, this yields the Hankel determinant B{N; —1/2). On the 
other hand, we have evaluated D{n,n — 1, N) in Lemma |1^. Then by dividing the results 
in ( [4.18 ) and ( l.l!j| ) by 2^(n — 1/2)^ and then setting n = 1/2, we obtain the expressions 
on the right-hand sides of ( [7.3|) and (|7.4| ). □ 

When solving the enumeration of rhombus tilings of a hexagon which contain the central 
rhombus, the determinant D{n,n -\- 1,N) was explicitly ||16[ (see the proof of Lemma 10 
and Proposition 14 in |jl6[) or implicitly evaluated (compare Proposition |ll]). If 

one adapts the preceding proof to this case, one obtains the following result. 
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Theorem 29. Let m be a positive integer. Then for 

B(N; 1/2) = det ((2^~'-^ - 

l<i,j<N ^ / T-jy 

there hold 



B{2m-1, 1/2) - ._,,,2,2 11 r9^^._ni 2^ 



24m-2 _ ly/ \^ 1 1 (2m + i - l)!y ^ (2m - 2z - 1) 

(7.5) 

and 



Proof. Consider the determinant D{n,n + 1, A^). Factor 2(n + 1/2) out of each column 
of D{n,n + l,iV), and then set n = —1/2. By the appropriate variant of ( |4.5|) and de 
I'Hospital's rule, this yields the Hankel determinant B{N; 1/2). On the other hand, by 
Proposition |Tl| (with replaced by + 1, / = n, and M = 2n — N) we know that 
D[n,n + 1, N) counts, up to a multiplicative constant, the number of rhombus tilings of a 
hexagon with sides A^ + l,2n — A^, iV + l,iV + l,2n — A^, A^ + 1, which contain the central 
rhombus. This enumeration problem was solved in |1^, 16|. If in the result we perform 
the according manipulations and then set n = —1/2, we obtain the expressions on the 
right-hand sides of ( |7.5| ) and ( |7.6| ). □ 



8. Concluding comments and open problems 

We conclude this article by pointing to open questions which are raised by this work. 

(1) In Corollary |^ it was demonstrated that, for M close to A^, the number of rhombus 
tilings of a hexagon with side lengths A^, M, N, N, M, N, which contain the rhombus above 
and next to the center of the hexagon, equals (|+r(A^, M)) T{N, M), where T{N, M) is the 
total number of rhombus tilings of the hexagon, and where r(A^, M) is a "closed form" ex- 
pression. (Magically, the value of 1/3 which appears here is, according to 0, Corollary 3], 
the exact proportion of the rhombus tilings that contain the central rhombus in the total 
number of rhombus tilings of a hexagon with side lengths 2n—l, 2n, 2n—l, 2n—l, 2n, 2n—l 
or with side lengths 2n, 2n — 1, 2n, 2n, 2n — 1, 2n.) As we mentioned in the Introduction, 
it is easy to derive many more such results, also for the central rhombus and the other 
two cases that were considered in Theorems Our proof, given in Section 2, consisted 
of starting with the expressions ( |1.2|) and ( p.. 3D and applying Zeilberger's algorithm to 
establish the simplification of the sum in these expressions when m and n are close. This 
is, unfortunately, not conceptual (as it just verifies, but does not derive the result), and 
therefore does not explain why these simplifications take place. The fact that apparently 
many more such results exist indicates that there must be a hypergeometric transforma- 
tion formula lurking in the background, which we were, however, unable to discover. (It 
is obvious that the sums in (|1.2|) - (|1.7|) can be written as very-well-poised yFg-series — 
see e.g. ( p.7| ) — and, by means of Whipple's transformation formula ( p.8| ), can therefore 
be transformed into balanced 4F3-series, to which, in turn, we could apply Sears' 4F3 
transformation formulas. However, it seems that this does not suffice to find the desired 
identity which would "explain" Corollary]^) 
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(2) Is it possible to find a uniform formula for the number of rhombus tilings of a hexagon 
with side lengths N, M, N, N, M, N, which contain an arbitrary (but fixed) rhombus on 
the "vertical" symmetry axis (i.e., the symmetry axis which runs in parallel to the sides 
of length M)? Recall that (as we mentioned already in the Introduction) in such a 
formula was found for the "horizontal" symmetry axis (i.e., the symmetry axis which cuts 
through the sides of length M). In contrast, here we encountered increasing difficulties in 
the proofs of our enumerations the farther we moved the rhombus which is contained in 
every tiling from the center. Recall that for solving our enumeration problems we needed 
to compute the determinants D{n,n — t,N) (see (|2.1| ) for definition) for t = 0, 1,2. For 
the case of a rhombus which is even farther away from the center, we would have to 
evaluate this determinant for even larger values of t. The increasing difficulties in doing 



this arise in Step 5 (compare the proof of Lemma |T2D of the computation. The previous 
steps. Steps 1-4, would prove that 



D{n, n — t,N) = (product of linear factors in n) ■ S{n, N, t), 



where S{n,N,t) is a polynomial of degree 2[(A^ + t)/2] (compare (|]2|), ( [4.20| ), and 
( [4.25| )). Thus, in order to determine S{n,N,t), the larger t becomes, the more eval- 



uations of S{n,N,t) at special values of n (or other informations about S{n,N,t)) we 
need. (The computations in |l^ have exactly the opposite behaviour: The farther the 



rhombus which is contained in every tiling is moved away from the center, the smaller 
in degree becomes the irreducible polynomial in the result.) Even worse, the larger t 
becomes, the more difficult it becomes to obtain these special values. (Remember, for 
example, the difficulty of evaluation of R{n, N) at n = 1 via Lemmas |25| and p6| .) 

That the problem that we considered here is at a different level of complexity than 
the problem in [|l^] is also indicated by the (partially conjectural) form of the asymptotic 
behaviour of the proportion of the rhombus tilings that contain this particular rhombus 
in the total number of rhombus tilings. While the asymptotic behaviour is totally smooth 
when the rhombus which is contained in every tiling is moved along the "horizontal" 
symmetry axis (see [|T^, Theorem 1.3]), the conjectured form |]TU|, Conjecture 1] of the 
asymptotics when the rhombus which is contained in every tiling is moved along the 
"vertical" symmetry axis behaves nonsmoothly. It is increasing for some time when the 
rhombus is moved away from the center, but at some point, when the rhombus enters the 
"arctic region" near the (top or bottom) corner, it becomes 1 and stays 1 from thereon. 
Thus, a formula for exact enumeration must, somehow, reflect this nonsmooth asymptotic 
behaviour. 

Is there a way to overcome these difficulties? 



(3) In Theorem ^ only c or d may be 0, but not a or b. In fact. Theorem ^ is wrong 
if a = or 6 = 0. But, apparently, not terribly wrong. Lemma ^ shows the evaluation 
of the determinant in ( |5.13| ) with a = b = 0, c = d = 2. Remarkably, the result is almost 
identical with the right-hand side in ( |5.13| ), the only difference being the polynomial in n 



of fourth degree in ( p.l|) . In fact, computer experiments suggest that a much more general 
result holds. 
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P{n;c, d), 

(8.1) 



Conjecture. For positive integers integers c, d there holds 

l<i,j<n 

1)! (c- 1)! (rf- 1)! {d-iy: 
{c + d-iy. 

c + i-l) (d + i-l) (d + i-l) {c + d + i~2) 
c + d + 2i-3) {c + d + 2i-2f {c + d + 2i-l) 

where P{n; c, d) is a certain polynomial in n of degree 2(c + d — 2). 
Furthermore, for positive integers integers b, c, d there holds 

det (5^+^-1 {B + 1)5-1 i-B + l),_i i-B + l)rf_i) 

l<i,j<n 

^ , ^^n(n-l)/2 / (c-1)! (& + C-1)! (rf-1)! (& + ^-l)! 

^ ^ V (6 + c + rf-l)! 

z (c + i-l) {b + c + i-l) (d + i-l) 
(b + c + d + 2i-3) (b + c + d + 2i-2)'^ 

(b + d + i-1) (b + c + d + i-2)Y~\ 
" (5 + + ^ + 2.-1) j ^(-'^'-'^)' (8-2) 

where R(n; b, c, c?) a certain rational function in n, which can he written with a numer- 
ator of degree c + d — 2 and a denominator of degree b — 1. 




In principle, our approach of proving Lemma [25| (the special case c = d = 2 of (|8TT| )), 
which consisted of using linearity of the determinant in order to break it into several 



pieces, to each of which we could either apply Theorem ^ or Theorem |2^, should make a 
proof of the above conjecture possible. However, serious difficulties have to be expected 
in actually doing the calculations, in particular, when working through a generalized form 
of Lemma ^ We believe that, in view of the simplicity of the result ( |6.1| ) and of the 
conjectured results ( ^.1| ) and ( |8.2| ), there must be a more elegant way to attack these 
Hankel determinant evaluations, in particular, if one also desires to obtain explicit forms 
for the polynomial P(n; c, d) and the rational function R(n; b, c, d). 

Note. Since first versions of this article were distributed, Use Fischer ("Enumeration of rhombus tilings of a 
hexagon which contain a fixed rhombus in the centre", preprint, nath/9906102 ) generalized Theorems ^ and ^ to 
arbitrary semiregular hexagons. 
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